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ABSTRACT 

Clines  ( nonuni  Form  spatial  1: str i but  ions  in  the  genetic  composition  of 
population  in  equilibrium'  ate  often  modeled  by  nonconstant  solutions 
u ! x )  c  [0,1]  of  (D(x)u')'  +  h(x,u)  =  0,  -»  <  x  <  *>,  where  h  satisfies 
h(x,0)  =  h(x,1)  =  0/  and  D  is  often  taken  to  be  identically  1.  The 
functions  n  and  h  have  interpretations  in  terms  of  mobility,  carryma 
capacity  and  natural  selection.  We  define  dines  a?  stable  solutions 
satisfying  u(-®)  =  0,  u(“)  =  1.  All  past  analyses  of  dines  have  consiiere 
the  case  when  (say)  0  is  the  favored  state  for  large  negative  x,  and  1 
for  laroe  positive  x  ( i  .<?. ,  ^  h(x,u)du  changes  sign  from  negative  to 

oositive  as  x  increases  from  -•»  to  ®) .  In  this  paper,  however,  we  assem 
that  the  state  0  is  favored  fm  all  x,  although  both  n  and  1  are 
stable  as  uniform  states.  A  number  or  conditions  are  aiven  which  ensure  th 
existence  of  stable  dines,  or  their  a*. dons  in  the  bounded  habitat  case. 
Conditions  are  also  gi”en  which  ensure  the  nonexistence  of  dines.  T'-e 
concent  of  stability  is  with  reference  to  the  corresponding  nonlinear 
diffusion  equation,  an!  is  use  1  in  a  special  technical  sense. 

AMS  (MOS)  Subject  Classifications :  R 2 A 1 d ,  D2A10,  34B15 ,  35K5R. 

Key  Words:  Clines,  stability,  population  aenetics,  ecology,  subsol ut i , 
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CLINES  INDUCED  BY  VARIABLE  SELECTION 
AND  MIGRATION 


P.  C.  Fife  and  L.  A.  Peletier 

1 .  Introduction. 

Clines  are  nonuniform  spatial  distributions  in  the  genetic  composition  cf 
a  population  in  equilibrium.  They  have  been  studied  a  great  deal  by 
population  geneticists  (see  for  example  [7,  9,  15,  16,  19]),  often  within  the 
context  of  the  simple  class  of  equations 

.2 

— s(x)f(u)  =  0  ,  (1.1) 

dx 

m  which  u  is  a  scalar  variable  representing  some  qene  frequency,  f  a 
selection  pressure  mechanism  tending  to  drive  the  population  to  certain 
preferred  states,  s  a  measure  of  the  selective  intensity,  and  the  second 
derivative  of  u  a  term  designed  to  account  for  spatial  migration. 

Outside  the  context  of  population  genetics,  equation  (1.1)  has  been  used 
in  models  designed  to  describe  certain  phenomena  in  ecology  [12]. 

We  suppose  that  f(0)=f(1)=0,  so  that  u  =  0  and  u  =  1  are 
solutions  of  (1.1).  These  solutions  represent  the  situations  where  the  gene 
in  question  is  either  entirely  absent  or  everywhere  present.  We  shall  say 
that  u  =  0  (u  =  1)  is  favored  at  the  location  x  if 

,1 

s(x)  J  f(u)du  <  n  (>  0) 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-004 1  and 
the  National  Science  Foundation  under  Grant  Nos.  MCS  78-04443  and  MCS  78- 
02158. 


One  reason  for  this  terninolouy  is  the  fact  that  in  the  case  when  s  is 
constant,  travel  lino  waves  of  the  corresponding  nonlinear  diffusion  eauation 

'it  -  uxx  +  s  f(u)  (1.2) 

always  move  in  a  direction  so  ♦•hat  at  ear  fixe’  x,  u(x,t)  tends  towards  the 
favored  state. 

Another  reason  for  this  terminolcrc.  is  that  it  is  consistent  with  the 
notion  of  relative  fitness  in  population  aenetics.  In  that  context,  the 
derivation  of  (1.1)  proceeds  from  a  s inale- 1 ocus ,  two  allele  model  (A  and 
a),  with  u  the  frequency  of  the  a-allele.  In  the  special  case  when  the 

fitnesses  of  the  genotypes  AA,  Aa  and  aa  are  independent  of  population 

density  and  genotype-f requencies ,  and  are  denoted  by  1  +  s^,  1,  1  +  Sj 
respectively  (s.,,  Sj  £  R) ,  sf  becomes 

sf(u)  =  ru(  1  -  u){(s1  +  s^lu  -  }  ,  r  >  0 

and  it  is  easi.lv  verified  that  the  state  u  =  f)  is  "favored"  in  the  ahove 

sense  if  and  only  if  s1  >  $2*  In  this  case,  travelling  waves  tend  to 
eliminate  allele  a,  and  drive  u  to  0. 

We  define  dines  to  he  stable  nonconstant  solutions  u(x)  of  (1,1).  x? 
t'e  domain  -  -  the  habitat  -  is  the  entire  real  line,  we  shall  for 
definiteness  require  that  u(-°°)  -  0  and  u(+<*>)  =  1.  r>n  the  other  hand,  if 
is  a  hounded  interval:  =  (-1.1),  we  reouire  Neumann  (uy  =  n) 

conditions  at  the  endpoints  and  u(1)  >  uf-1).  my,P  notion  of  stability  is 
here  to  he  understood  with  reaard  to  equation  (1.2).  The  precise  definition 
of  stability  to  he  used  is  critical,  and  we  shall  explain  it  at  the  end  of 


this  section.  The  reason  for  including  stability  in  the  definition  of  a  cline 
stems  from  the  fact  that  unstable  solutions  of  (1.1)  are  not  likely  to  reflect 
anything  seen  in  the  natural  world. 

Past  analyses  of  dines  [3,  6-9,  15,  16,  19]  have  all  been  devoted  to  the 
case  when  u  =  0  is  favored  for  large  negative  x  and  u  =  1  is  favored  for 
large  positive  x.  In  this  paper  we  shall  suppose  that  u  =  0  is  favored 
everywhere  in  the  habitat: 

s  >  0  in  R  and  /J  f(u)du  <  0  , 

and  we  ask  whether  it  is  still  possible  for  dines  to  exist.  If  s  is 
constant,  the  answer  is  no.  For  a  bounded  habitat  this  follows  from  [1],  and 
if  r.  =  R  it  follows  from  results  in  [10]  and  elsewhere.  However,  we  shall 
show  that  provided 

(i)  s(x)  is  allowed  to  vary  in  certain  ways,  and 
(ii)  the  selection  mechanism  is  bistable,  i.e. 

f ’ (0)  <0  and  f ’ ( 1  )  <  0  , 

a  dine  can  indeed  exist.  In  terms  of  the  example  from  population  aenetics, 
bistability  means  we  are  in  the  underdominant  case,  i.e.,  s1  >  0  and 

s  p  >  0  . 

Ana  in  in  the  nondation  nenoiic  context,  (1.1)  is  derived  under  the 
ars'j'rtions  that  the  total  population  density  is  constant  in  space  and  time, 
that  *  he  migration  rate  is  likewise  constant,  that  no  drift  (i.e.,  no 
preferred  direction  of  mijratior)  j  present  in  the  latter,  and  that  the 
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re la Vive  fitnesses  of  the  genotypes  are  constant.  These  restrictions, 
however,  can  be  dropped  and  (1.1)  replaced  by  a  more  general  equation.  These 
centralizations  are  examined  in  this  paper.  For  example  in  sections  6  and  7, 
ve  extend  our  analysis  to  allow  for  variable  total  density  and  migration  rate 
'but  we  continue  to  assume,  for  simplicity,  that  no  drift  is  present),  and 
find  that  spatial  variation  in  the  carrying  capacity  of  the  environment  and  in 
the  migration  rate  may  also  be  sufficient  for  the  existence  of  a  cline. 

The  biological  picture  we  shall  have  in  mind  is  the  following.  The 
individuals  of  the  three  genotypes  AA,  Aa ,  aa  operate  under  the  same 
migration  rules:  more  specifically,  there  is  an  infinitesimal  variance 
(mobilitv)  V(x)  in  the  migration  rate,  common  to  all  three  types,  and  in  all 
cases  the  drift  is  zero.  (The  effects  of  genotype-dependent  migration  were 
studied  in  [14].]  The  carrying  capacity  k  of  the  environment  is  assume  to 
be,  to  first  order,  independent  of  the  population’s  genetic  composition,  and 
is  a  given  function  k(x).  When  the  relative  fitnesses  of  the  genotypes  vary 
in  space,  the  selection  term  s(x)f(u)  should  be  replaced  by  some  function 
hfx,u).  Under  all  these  conditions,  the  appropriate  generalization  of  (1.2) 
is  [5,  p.M] 


JU 

3t 


V  (  x )  V  (x) 


yy  (V" ( x)k" ( x)  yy)  +  h(x,u) 


(1.3) 


Other  mi  oration  rules,  such  as  qiven  in  (5,  (1.20)]  could  be  handled  hv 
our  technique;  but  the  results  obtained  here  suffice  for  the  purpose  of 
i 1 lustrating  the  effects  of  variable  mioration  and  carryina  capacity. 

We  shall  oive  various  sufficient  conditions  for  the  existence  of  dines; 
typically  they  involve  the  statement  that  the  function  s,  V,  or  k 
experiences  a  sufficient  drop  over  some  finite  zone  in  the  habitat.  However, 
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not  only  does  the  magnitude  of  the  drop  in  this  zone  have  to  be  large  enough; 
but  also  the  gradient  of  the  function  in  question  has  to  be  steep  enough  over 
part  of  this  zr  :  The  function  should  remain  relatively  large  over  that  part 

of  the  habitat  lying  to  the  right  of  the  zone,  whereas  its  values  to  the  left 
are  not  so  important.  The  intuitive  role  of  this  zone  is  that  it  separates 
two  communities,  which  are  more  or  less  in  states  u  =  0  and  u  =  1 
respectively.  In  the  case  when  the  function  experiencing  the  drop  is  the 
mobility  V,  the  zone  would  clearly  act  as  a  barrier  tending  to  isolate  the 
two  communities ,  and  so  providing  for  their  coexistence  in  different  states. 
Our  results  show  that  a  similar  variation  in  carrying  capacity  or  strength  of 
selection  can  function  eoually  well  as  a  barrier.  Matano  [13]  has  shown  that 
in  higher  dimensions,  the  shape  of  the  habitat  ft  is  relevant  to  the 
existence  of  dines,  domains  with  narrow  middle  sections  favoring  their 
existence  through  a  type  of  barrier  action. 

We  shall  give  a  brief  description  of  the  results  of  the  paper.  The 
function  u  will  assume  values  only  in  the  range  0  i  u  <  1.  The  following 
hypotheses  on  f  and  s  will  be  made  throughout. 

H  :  f  e  C2([0,1]);  f(0)  =  f(a)  =  f(1)  =  0  for  some 

a  e  (0,1);  f'(0)<0,  f'(1)<0;  f(u)  <  0  on  (0,a)  , 

f(u)  >  0  on  ( a ,  1  ) ;  and  f(u)du  <  0 

Remark :  This  last  inequality  means  (since  s  >  0)  that  the  state  u  =  0  is 

favored  for  all  x. 

H  :  s  is  piecewise  continuously  differentiable  on  il  where  i.  =  R  or 

[-1,11; 

i  n  f  { s  (  x )  :  x  <:  ,1 )  -  s  J)  >  0 
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Theorem  5  and  Corollary  5.1  give  conditions  on  h,  in  the  form  of 
inequalities,  which  imply  the  existence  of  a  cline.  The  principal  require¬ 
ments  are  that  h(x,u)  >  s(x)f(u),  where  Hf  and  Hs  are  satisfied,  s  is 
sufficiently  small  on  some  interval  (0,B),  s  >  1  on  (A,°°)  (where  A  >  B) , 
s  S  1  on  ( B , A) ,  B  and  A  are  sufficiently  large,  and  A-B  is  sufficiently 
small.  The  internretation  in  that  (0,R)  is  a  "barrier”  region  between  tbe 
two  "communities"  ana  (A,00),  and  ( B ,  A)  is  a  "transition"  zone. 

Finally,  in  Theorem  6,  it  is  shown  that  s'  >  0  implies  that  the  cline 
is  strictly  increasing. 

Results  for  a  bounded  habitat  and  for  variable  V  or  k:  An  existence 
condition  analogous  to  Theorem  4,  but  for  a  bounded  habitat,  is  given  in 
Theorem  7.  Finally,  conditions  for  nonexistence  and  for  existence  of  dines 
are  /iivan,  for  cases  in  which  V(x)  and  k(x)  are  variable,  in  Theorems  8  - 
ig.  ?'■  e  conditions  are  similar  to  the  conditions  on  s  in  the  foreooincj 
r/~-,  lit/  . 
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Results  o*  the  tyr<“  presented  i  n  *"bts  paper  were  conjectured  to  he  true 
>nn  (personal  communication  and  M2’).  In  M2],  the  int  erpretaticr 
wit-hip  the  context  of  the  spatial  variation  in  composition  of  ecolooica 
ties,  rather  than  the  context  described  here, 
have  assumed  throughout  that  ;  ^  f(u)du  <  0  (or  the  analoaous 
r<  :r.  on  h),  so  that  0  is  the  favored  state  for  all  x.  As  rent  i  one  : 

h-*f  •  ,  this  is  a  major  distinction  from  previous  treatments  of  dine*’ ,  v,k»r 

♦  he  f a vo rod  state  was  assumed  to  change  within  the  habitat.  It  aces  wi tho.jr 
cv  i'll  that  if  j  f  du  =  0,  then  a  cline  exists  even  in  a  homogeneous 
habitat;  hut  it  is  structurally  unstable,  since  a  slight  change  in  the 
■unction  f  will  destroy  it.  To  ensure  that  a  cline  exists  when  _  e  cm  - 
and  that  a  cline  continues  to  exist  when  f  is  perturbed  by  any  small  amour, 
it  suffices  to  have  s,  V,  and/or  V  vary  in  ways  described  by  our  tb.eorv: 
we  shall  not  pursue  the  details.  The  idea  is  that  existence  conditions  ray 
given  which  remain  valid  when  these  fun.-Mc's  are  altered  siinhtlv  but 
a  rb  l  trar  ily . 

Penally,  let  us  explain  the  notion  of  <-t ability  which  we  s“all  use 

■  :  ,;v;hout  the  paper.  It  is  designer,  so  that  th»*  following  statement  is  ; 

Proposition .  Let  £  be  a  station, irv  su!  sol  u*- ;  m.  ard  .*  a  Ptationarv 
uupersolution  of  (1.3)  together  with  appropriate  boundary  conditions,  such 
;  p  {  '  »  ,  and  neither  is  an  exact,  solution.  Tv  er.  there  exists  a  staMi 

■  r  r  i  tin  solution  -A  of  (1.3)  between  £  an  1 

is  proposition  is  the  basis  Fr  r  ri.ur)v  all  *  hr  ov*,s*ense 
*•  ■  ■  •  ■:-,t *  .  A  stationary  subso  a*-  s  o-  .  of  '  1  .  '  .  ■  i  <  r  i  •'  •  •  1  he  a 

-  *  *  i  ■  . , n u s  f  urict  i  on  on  ,  *■  ucv'  1 1  a  *  ‘  so -re  r  ..... .  i •  ■  ' ■  r  ;  *  "v"  '  *  ■ 

is  t  w  i  re  nu!  ‘ :  •'  ‘  .  r.  *•  . . .  ■ 

~C  i  er 


c  it  ■ 


nM  >  o 


N  [  u]  being  the  expression  on  the  right  of  (1.3),  whilst  at  each  endpoint  •" 
in  the  interior  of  2 ,  £'(?  +  0)  >  £'(C  -  0).  For  2  bounded,  we  also 
require  £'  <  0  at  its  left  endpoint,  the  opposite  inequality  holding  on  the 
right.  An  analogous  definition  holds  for  supersolutions. 

If  the  habitat  J2  is  bounded,  we  shall  define  stability  in  the  usual 

O 

C  sense.  It  was  shown  by  Matano  [13]  that  the  proposition  is  true  when 
stability  is  interpreted  this  way.  If  the  habitat  is  unbounded,  the 
proposition  has  only  been  proved  for  the  following  weaker  notion  of  stability 
[2]  . 

Let  y :  R  +  (0,1)  be  a  solution  of  (1.1),  and  let  T  be  a  family  of 
functions  (defined  below)  with  prescribed  behavior  for  large  ! x I . 

Definition.  ^  is  said  to  be  stable  if,  for  every  ^  e  T  with 

<  y  <  y^,  there  exist  elements  £  ^  ,  y'  ^  e  T  with  ^  <  $  <  y' ^  such  that 

if  y ^  <  u^  <  y 2  then  the  solution  u(t;u0)  of  (1.1)  which  satisfies 

u(0;uq)  =  u0  has  the  property 

lim  inf  u(t,u^)  <  lim  sup  utt.u^)  <  y'2 

£-►00  £>CD 

Here  all  the  inequalities  are  understood  in  the  pointwise  sense. 

To  complete  our  definition  of  stability,  we  need  only  say  which  function 
belong  to  the  family  T. 

Definition:  A  function  ^  e  C1 (R;  [0,1] )  is  said  to  belona  to  T  if 

(i)  y ( x )  is  an  exact  solution  of  (1.3)  for  laroe  Ixl ,  and 

(ii)  , ( -* )  =  n ,  v ( °° )  =  1 . 
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Within  our  level  of  general ity ,  it  will  not  necessarily  follow  that  the 


stable  solution  we  obtain  is  asymptotically  stable,  in  that  small 
perturbations  of  it  return  to  it  as  t  *  °°.  In  fact,  in  exceptional  cases 
there  may  be  many  dines,  and  all  we  sav  is  that,  when  one  is  perturbed  by  a 
small  amount,  the  subsequent  evolution  of  u(x,t)  never  leaves  a  neighborhood 
of  the  cline  that  was  perturbed.  A  different  cline  in  that  same  neighborhood 
may  be  approached.  We  are  guaranteed,  however,  that  u(x,t)  remains 
between  £  and  ~  for  all  fine,  and  so  despite  the  possibility  of  many 
successive  perturba*  ons  through  the  course  of  time,  it  will  never  evolve  to  a 
spatially  constant  distribution. 

Clines  can  also  be  studied  from  the  point  of  view  of  discrete,  rather 
*-'~a r.  continuous,  selection-miqration  models.  A  result  of  Karlin  and  MacGregor 
'IV,  for  example,  implies  that  when  the  migration  rate  between  colonies  is 
snail  ar.i  each  colony  in  isolation  has  two  stable  equilibrium  states,  then 
-anv  dines  exist.  This  is  conceptually  in  accord  with  our  results,  which 
say,  very  rouohly,  that  certain  spatial  heterogeneities  favor  the  occurence  of 
dines.  Colonies  with  small  interchange,  when  conceived  in  a  geographical 
setting,  form  an  inherently  heterogeneous  model. 

The  authors  be"of it*V  greatly  from  discussions  with  T.  Nagylahi  and 
S.  Levin. 


Nonexistence  Theorems. 

We  consider  the  problem 


(I) 


u"  +  Xs(x)f(u)  -  0 
u(-“)  =  0,  u(+°°)  = 


x  £  R,  X  >  0 , 


(2.1) 


where  f  and  s  satisfy  the  hypotheses  and  Hg  listed  in  the 

Introduction . 

It  is  well  known  that  if  si  constant,  Problem  I  has  no  solution  in 
view  of  the  fact  that  f(u)du  /  0.  Thus  it  is  the  variation  of  s  which 

may  give  rise  to  the  occurence  of  dines.  In  the  following  theorem  we  give  a 
lower  bound  for  the  variation  of  s  which  is  necessary  to  sustain  a  cline. 

We  define 


si  =  inf  six)  ,  s ^  =  sup  s(x) 

R  R 


Kid  we  assume  that  s^  >  0. 
Theorem  1 ,  Let 


s 

s 


2 


> 


f  (  u ) du 
'  a _ 

f  (  u ) du 


(2.2) 


Then  Problem  I  has  no  solution. 

o-  .  f.  Suppose ,  to  the  '-ontrary ,  that  Problem  I  has  a  solution  <  . 
d  i:  ir’y  .*  ’  ( x )  >  0  for  large  negative  values  of  x.  Set 


xn  --  sap;  x  :  v'  1  (  f  )  >  0  on  (  -°°  ,  x )  } 


or  xfl  <  ,  in  which  case 


-1  1  - 


Multiplying  (2.1)  by  u'  and  integrating  over  (~“»xq)  we  obtain 


^(x  ) 

/  s(x(sf  )  )f  (sff  )  ds9  =  0  , 


(2.3  ) 


where  x(^)  denotes  the  inverse  function  of  ^(x).  Since  ^  '  >  0  on 
(-,x|  this  function  is  well  defined.  Since  t  ( x^ )  £  (a,1)  (2.3)  implies 

that 


/p  s(x(*  ))f(V')ds;  >  P  . 

Thus,  using  the  bounds  for  s  we  obtain 

s1  /*  f  (  v ) d g  +  s2  f(y')dy-  >  0 


or 


J1  f(g)d  y' 

f'v)d > 


(2.4) 


which  contradicts  our  assumption  about  Si/S2* 

Theorem  1  is  optimal  in  the  following  sense:  If  positive  numbers  s1 
and  s2  are  given  which  violate  (2.2),  there  exists  a  function 
s  :  R  +  [8l,8  1  such  that  Problem  I  has  a  solution.  We  shall  demonstrate 
this  in  the  following  proposition:  (see  also  [17]). 
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Proposition  1 .  Let  s 


R  ■*  R+  be  defined  by 


s(  x) 


fs  i  *  <  0 
\  *2  X  >  ° 


where  0  <  s1  <  s2  and  s^/s2  satisfies  (2.4).  Then  Problem  I  has  a 
solution . 


Proof.  We  construct  a  supersolution  u  and  a  subsolution  u  such  that 
u  4  u.  Then  the  existence  of  a  solution  y  such  that  u  <  ■f  <  u  is  ensured 
by  the  standard  theory  [18]. 

We  begin  with  the  construction  of  a  monotone  subsolution.  Let  V1  be 
the  solution  of  the  problem 


J/"u"  +  s1f(u)=0  -“  <  x  <  0 

\  u( -«)  =  0,  u( 0 )  =  a 


(2.5) 


and  let  V2  be  the  monotone  solution  of  the  problem 

(  u"  +  s-,  f(u)  =0  0  <  x  <  “ 

/ 

l^u  ( 0 )  =  a ,  u  ( ” )  =  1 


(2.6) 


Both  solutions  exist,  and  can  be  constructed  implicitly  by  quadrature.  In 
fact,  multiplying  (2.5)  by  u'  and  integrating,  we  obtain 


i  {u' 


(  x)  }  +  s 


/If  f(s)ds 


1  J0 


0 


which  qi ves  u'  in  terms  of  u  for  u  e  (0,a] ,  since  the  integral  is 
negative  in  that  range.  Define 
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Then  u  will  be  a 


subsoltition  of  Problem  I  if 


—  + 

u* (0  )  «  u' (0  ) 


It  follows  from  (2.51  that 


{V  1  ’ i 0  >  J7  =  -s1  jp  f ( v)dv 


and  from  (2.6)  that 


1V2’(0)}2  =  +s2  f(v)dv  . 


Thus ,  by  (2.4) 


{V  1 (0)}“  «  {V  ' (0)}' 


which  establishes  (2.7). 

For  a  supersolution  we  take,  for  some  5  £  the  Action 


u(  x) 


; 


w(  x) 


X  >  f, 

x  <  £ 


where  w  is  the  solution  of  the  problem 


This  problem  has  a  unique  monotone  solution  (use  the  method  outlined  above 
and  Hjr) .  By  choosing  £  <  0,  we  make  u  a  supersolution  of  Problem  1,  and 
by  choosing  (-£)  large  enouqh  we  can  ensure  that  u  >  u.  This  completes  the 
proof . 

In  the  next  theorem  we  limit  the  class  of  functions  s  for  which  we  can 
expect  dines  even  further. 

Theorem  2.  Let  s  be  a  nonincreasing  function  of  x.  Then  Problem  I  has  no 
solution . 

Proof .  Suppose  again  that  Problem  I  has  a  solution  q ,  and  define  xQ  as  in 
the  proof  of  Theorem  1.  Then  we  obtain  (2.3)  which  we  now  write  as 

£ 

-/*  s(x(v'))f(;)d;=  /  0  s(x(y))f(vOdq  (2.8) 

with  y  q  =  v(  Xq  ) .  Since  v  is  monotone  on  (-®,xo)  and  >?  (  Xq  )  >  a,  there 

exists  a  unique  point  £  e  )  such  that  c(£)  =  a.  Now  utilizing  the 

monotonicity  of  s  we  deduce  from  (2.8)  that 

y 

-s(£)  /*  f ( v ) dc  <  s(?)  |an  f(v )dc 


and  hence 


-1 

'  n 


f  (  c )  d  v- 


f(y)d  y-  >  0  . 


ntradi ots  our  assumption  about  the  si  an  of 


,1 

1  0 


r 


-If- 


Thus  for  there  to  exist  a  cline,  it  is  necessary  that  s  be  increasina 
at  some  point  of  the  domain.  In  fact,  as  we  shall  see  as  a  corollary  of  the 
following  more  general  nonexistence  theorem,  this  is  still  not  enough:  it 
needs  to  increase  sufficiently  fast. 

To  prove  this  we  consider  the  problem 


(II) 


i  u"  +  h(ex,u)  *  0 
)  u(-«)  =  0,  u ( 00 ) 


1 


where  we  make  the  following  assumptions  about  h. 

Al.  h  C  C1 (*  x  (0,11),  h(x, • )  c  C2([0,1)) 

A2 .  | h  |  t  1 ,  | h  1  t  M  ,  | h  |  <  M  uniformly  on  R  x  [0,1]. 

x  u  1  uu  2 

A3.  h(x,0)  =  h(x,a(x))  =  h(x,1)  =  0  for  all  x  e  R,  a  e  C(R), 

0  <  a(x)  <  1  for  all  x  c  R;  h(x,«)  <0  on  (0,a(x))  and 
h( x  ,  • )  >  0  on  ( a( x) , 1 ) . 


A4. 

A5. 


h  (x,0)  t  -a  <  0,  h  (x,1 )  <  -a  <  0  for  all  x  t  R. 
u  u 

h(x,u)du  <  -3p  <  0  for  all  x  e  R. 


Note  that  Problem  II  is  equivalent  with  the  problem 


f  u"  +  A  h(x,u)  =  0  X  =  1 /e2 

(II')  (2.9) 

,  u(-°°)  =  0,  u(®)  *  1 


Theorem  3.  Let  h  satisfy  the  assumptions  A1-5.  Then  there  exists  an 

*  *  *  * 

c  >0  (X  >0)  such  that  if  etc  (X  >  X  )  Problem  (II)  (II')  has  no 

solution . 


Proof .  Suppose  «f(x,X)  is  a  solution  of  Problem  TI',  and  let  be  the 

smallest  root  of  the  equation 


-If- 


■£(x,X)  =  a(x) 

In  view  of  the  boundary  conditions,  and  the  assumptions  on  a,  such  a  root 
will  always  exist.  Note  that  £  will  depend  on  X. 

Before  we  proceed  further,  we  introduce  some  notation.  Choose  d  >  A 
such  that 


/ .  h(x,u)-du  <, 


-2p 


for  all  x  e  R 


and  define 


R(C)  =  {(x,u)  :  |x-£|  <  P,  5  <  u  <  1} 

To  begin  with  we  shall  show  that  there  exists  a  X  >  0  such  that  i r 
a  >  the  graph  of  <f  can  only  enter  R  through  the  bottom,  i.e. 

y ( x ,  X )  <  {  on  (-«,$-p]  .  (2. 10' 

Suppose  to  the  contrary  that  v'  enters  R  at  a  point  (£-p,u  ),  where 
u1  e  (5,1).  Then,  because  y'"  >  0  on  (£-p,£), 

v’(C-p)  <  ( 1  —  6 ) /p 


and  hence 


j'lZP  I  h  ( x  ,v*  ( x ,  X  )  )  |  dx  <  (1-5)  '(o\)  . 


(0.11) 
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1 

f 

i 


1 

I 


Since 


1 4—  h(x,v(x,A))|  *  I  h  (x,y(x,\))|  +  |  h  (x,^(x,A))|  |s?'  (  x ,  X  )  | 
dx  x  u 

s  1  +  M  (1-5)/p 

the  inequality  (2.11)  implies  that  at  each  x  e  (-“,£-p], 

h(x,^(x,A))  *  0  as  A  *  °° 


and  in  particular 


h(£-p,  y'(£-p,A))  -*■  0  and  A  ■+  °°  . 


Thus,  for  A  large  enough, 


u  1  =  v  (  £  -p  ,  A  )  <  5 


and  we  obtain  a  contradiction. 

We  shall  now  show  that  there  exists  a  A^  >  0  such  that  if  A  >  Aj ,  the 
graph  of  c  car.  only  leave  R  through  the  top.  This  contradicts  the 
boundary  condition  at  +°°  and  thereby  proves  the  theorem. 

We  multiply  (2.9)  by  2v’  and  integrate  over  (-“>,£).  This  yields 

i  ^ '  ( £  ) } 2  =  2  A  ^(;)lh(x(;),<)M;>  2A  |h(x(*)  ,*)  |d  *  » 


-in- 


,,pre  is  the  inverse  of  <x).  Let  X^.  Then  |x(y'l- 

f*  aT  and  it  follows  that 


51  <  P 


on 


v  . .  {  r  ,  > 2  >  2  A  ;fC)|h(C,v'lMy' ~  2Apa(5)suplhxlsup|x(^)  -  Cl 


>  2A  r,(^  |h(C.  v)l<3y'  -2Xpa(C)  . 

"  0 


(^,C+p1  <  we  now  aet <  as  ^-onC!  as  ^  >  ^ 


UMx)}2  =  iC(C);2  -  2A  ra  h(  x(  v ) ,  s?)  <V 


>  -2 A  Jr  MC.vldv 


;  -2 Apa  -  2A  /fX)  h(5.v">d< 


2Xp ( 1 -a ) 


>  +  2 X { -/ r  h(C,v)dr  -  pl 


>  2a ( 2p-p ) 
=  2Ap  . 


Thus,  choosing  A.,  =  1/2p52,  we  achieve  that 


(x)  >  1/5  x  >  5 


vl  hence,  that  the  graph  of  leaves  P  through  the  top. 


cenar 


arv.  Theorem  3  mav  he  regarded  as  a  partial  converse  to  Theorem  fi.1  of 
,n  wMch  thc  existence  of  a  transition  layer  solution  is  established  for 


*-t.  r-  *riia  *■  i  on 


+  hfx.u)  ~  h 


x  ‘  R 


-r- 


in  which  h(x,0)  =h(x,1)  =  0,  h  (x,0)  <  -<  h  (x,1)  <  -<  for  some  <  >  0 

u  u 

and  in  which  it  is  assumed  that  for  some  xQ  £  R, 

/p  h(  xQ  ,u)  du  =  0  . 

In  our  terminology,  a  transition  layer  solution  is  a  special  type  of  cline, 
one  on  which  the  change  is  primarily  concentrated  in  a  small  interval.  It  is 
the  last  condition  which  is  violated  here,  and  indeed  we  find  that  no  cline 
exists  for  small  enough  £. 

We  now  return  to  our  original  problem.  Consider 

(  u"  +  s( ex) f ( u)  =  0 

U  ) 

'  l  u<-“)  =  0,  u(“)  =  1  , 

where  s  and  f  satisfy,  respectively,  Hg  and  Hf.  In  addition  we  assume 
that  s  £  C1 (R) ,  and 


ls'(x)t  <  m,  xeR 


*  * 

Corollary  2.  There  exists  an  e  >  0  such  that  if  £  <  e  ,  Problem  I  has 
no  solution. 
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^y 


i 


*l:L 


3 .  Existence  Theorems. 

Having  established  a  number  of  necessary  conditions  on  s  for  there  to 
exist  a  cline,  •>  shall  now  give  two  sets  of  sufficient  conditions. 

In  the  first  theorem,  we  describe  the  effect  of  a  local -'zed  inhomogeneity 
in  an  otherwise  uniform  habitat.  Thus,  we  assume  that  s  is  given  by 


s(x,u;£)  =  {l  +  u9(x,£)}  1  u  >  0,  £  >  0  , 


(3.1) 


where  0  z  C  (R)  has  the  following  properties: 

HI  0(x,C)  =  0  for  I  x |  >  £ 

H02  0 (x,£ )  >  0  for  |xl  <  £  (=0) 

H.3  Max{0(x)  :  |x|  <  £}  =  1 

O 

Thus,  2£  denotes  the  width  of  the  inhomogeneity,  and  u  is  a  measure  of  its 
"strength" . 


Theorem  4.  Let  f  satisfy  Hf  and  let  s  be  of  the  form  (3.1)  in  which 
satisfies  the  assumptions  H  .  Then,  if 


.  _  2  a  .  f  t  *  *•  1 

Xt,  >  ~z~  {  f(r)dr}  , 

8  ‘  a 


(3.2) 


there  exists  a  number  u  such  that  if  u  >  M  ,  Problem  I  has  a  solutior 

* 

U  depends  on  X  and  £  and  tends  to  infinity  if  either  X  •+  <*>  or 
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Remark:  We  have  normalized  the  function  sfx)  so  that  its  maximum  is  equal 


to  one.  This  can  be  done  by  properly  scaling  the  variable  x.  By  (3.1),  we 
are  assuming  that  s  is  identically  one  for  |x|  >  £  but  experiences  a  drop 
in  the  center  zone  |x|  <  .  Tt  will  be  clear,  however,  that  the  proof  works 
as  well  if  s  is  any  other  positive  constant  for  x  <  Thus,  let 

be  a  smooth  monotone  function  with  x(*)  H  a  >  0  for  x  <  -£  and  =  1 

for  x  >  Z, .  Multiply  the  function  on  the  right  of  (3.1)  by  x(y^  and  use  the 
resulting  expression  for  s.  Then  Theorem  4  and  its  proof  are  still  valid. 

Proof  of  Theorem.  We  construct  a  sub-  and  a  supersolution.  The  subsolution 
v  we  are  going  to  construct  will  consist  of  three  pieces:  vj ,  v2,  and  v3 . 
For  v^  we  take  the  solution  of  (2.1)  on  (£,,“)  which  satisfies  the  boundary 
conditions  v  ( £)  =  a  and  v3(00>  =  1.  This  solution  clearly  exists,  and  it 
is  unique.  For  v1  we  choose  v1(x)  =  P  on  (-°°,p],  where  pc  (-f,,C)  is 
still  to  be  determined.  Observe  that 

{ v  '  f r  )  }  2  =  21  '1  f  (  r)  dr  . 

3  •  .3 

Let  c  ( u)  >  or,  [0 ,  a]  .  Then  we  choose  for  Vo  the  solution  of  the 

problem 

;  U"  -  1st  X,  U,  ",  )M  on  (p,S) 

(3.3) 

•-(')  -  a,  u'  '  %)  =  v,‘  (D  . 

Then  vp  is  1  subsolution  of  (2.11  on  ( p  ,•] )  •  Wo  shall  show  fhat  if  (3.2)  is 
sat  isf  iei  and  large  enough,  vV  >  0  on  (n,y)  and  v^Cp)  =  P  • 


Cloarly,  since  u"  >  0  by  (3.3)  a  necessary  condition  is 


v3'(rJ  >  a/ (20 


which  implies  (3.31.  Moreover,  integrating  (3.3)  we  find 

ft"  r}  v  f 

u'(l)  -  u'(x1  <  X«  -  , - —  =  £  .  (3.4) 

'  -i  1  +  uc ( x  ,C  ) 

Thus  integrating  once  again,  we  obtain 


a  >  {v3  '  (;>  -  e}  (f,-p) 


whence  p  *  -  t,  if 


This  can  clearly  he  achieved  by  choosing  u  >  U*  for  some  large  u*.  It  is 

★ 

clear  however  fro*~  P.4)  that  as  X  “  or  £  ♦  »,  u  *  »  as  well. 

The  super solut ion  w  is  constructed  out  of  two  pieces  w1  and  W2« 

Here  w1  is  the  solution  of  (2.1)  on  (-“,-£),  which  satisfies  the  boundary 
conditions  w  (-*>)  -  0,  w  ( —  1 )  =  1.  In  view  of  the  integral  condition  on  f 
this  solution  exists,  and  is  unione.  por  w^  we  choose  w^(x)  =1  on 
/  1  • 

construction  v<  x)  <  w  ( x)  on  R.  Hence  there  exists  a  solution  u 
of  Prnhle"-  t  such  that  >•  <  u  w. 


-.'3- 


To  construct  other  sufficient  conditions  for  the  existence  of  a  stable 


cline,  we  envisage  a  pair  of  piecewise  continuous  functions  H(u)  and  H(u) 
satisfying 

H(0)  =  H( 1 )  =  0  ,  (3.5) 

/k  H(u)du  <  0  for  all  b  e  (0,1]  ,  (3.6a) 

1  H(u)du  >  0  for  all  r>  £  [0,1)  ,  (3.6h) 

b  — 

together  with  solutions  U(x),  U(x)  of 

V"  +  H(U)  =  n,  U(®)  =  1  ,  (3.7) 

U"  +  H(U)  =  0,  TJ  ( -«> )  =  0  .  (3.3) 

Representing  the  solutions  in  the  usual  way,  one  can  see  that  U  must  vanish 
at  a  finite  value  x,  which  may  be  chosen  arbitrarily,  and  U  -  1  at  a 

finite  value  x.  We  consider  U  to  be  defined  only  on  [x,co)  and  U  only 

on  (  -00 ,  x  ]  . 

Theorem  5.  Let  H,  H,  TJ,  U  be  as  ahove,  with  U  (  x )  4  U(x)  where  they 
are  both  defined.  Assume  h(x,u)  satisfies  and 

h(x,tT(x)  )  >H(n(x)>  for  Uc  [0,11  ,  (3.na) 

h(x,U(x))  5  H(tT(  X)  )  ,  for  n  c  [0,1]  .  (  ' . n  b ' 


Then  there  exists  a  stable  solution  of  Problem  (II)  with  z  ~  1  . 

Proof :  The  function  u(x)  =  Max{0,U(x)j  is  a  subsolutir.r. ,  and 

u(  x)  =  Min{1,-.  ’’  ;  3  supersolution.  The  conclusion  follows  immediately. 

The  followina  coronary  treats  the  case  when  h  can  be  compared,  in  a 
certain  way,  with  functions  of  the  form  s(x)f(u),  where  s  and  f 
satisfy  H  and  Kr  ,  r«*srx»ct  !  ye  1  v .  We  picture  the  habitat  as  heine  divided 
into  four  regions,  the  Mrst  two  hsino  (i)  the  interval  (~a3,0),  where  the 
population  exists  in  a  s  ate  near  u  -  0;  and  (ii)  the  interval  ( A,® ’•  for 

some  A  '  0,  where  u  is  near  1.  The  interval  (0,A)  between  these  two 
population  states  is  divide.)  into  (iii)  a  "barrier"  region  (h,E),  wherein 
the  selection  strenoth  is  rrcuired  to  be  small  enough,  and  ( iv)  a  "transition 
interval  ( B , A)  of  length  T  =  A-B ,  which  is  simply  assumed  to  be  short 
enouah.  In  region  (ii),  the  selection  strength  must  be  larqe  enouah  but  no 
strona  reauirements  are  imposed  upon  it  in  (i)  and  (iv). 

The  critical  ratio 


V 


' 1  f ( u ) du 
•  a 

-/*f (u)du 


(3.10) 


will  play  a  crucial  role. 

Corollary  5.1.  Let 
f  satisfy  Hs  and  H ^ , 
numbers  B,  T,  and  o. 


Mote  that  0  <  y  <  1  . 

r  <  y.  Suppose  h(x,u)  f  s(x)f(u)  where  s  and 
respectively,  and  in  addition  for  some  nonneaat: 


S.(x)  *  C,  X  (-“,0) 

s(x'  •  ,  XS  ( 0 , B) 


s ( x )  *  1 ,  x  ( P , P+T 1 


(  3  .  1  la-.!' 


i 


six)  ■  1  , 


x  c  f  B  +T  , 00  ) 


Also  for  some  neaative  number  -X,  we  assume  that 

h(x,u)Ys(x)f(x)  ,  X  <  -X  , 

where  e  and  f  also  satisfy  !!,_  and  II,-,  0  <  s1  *  s(x)  3  for  x  <  -X , 

and 

15  1 

—  >  •-  ,  (3.12) 

s 

Y  beinq  defined  by  ( 3 . 1  ('  >  with  *  replaced  by  f.  Then  there  are  numbers 

*  * 

B  (y  ,0) ,  T  <y,c)  such  that  if 

-1/2*  - 1  /->  *  -1/1*  * 

B  >  X  /ZB  ,  T  A  '  T  ,  B  +  T  >  A  ~(B  +  T  )  ,  (3.13) 

there  exists  a  stable  solution  of  Problem  (II'). 

Pemark .  Formulas  for  suitable  F  and  T  will  be  qiven  in  Proposition 
5.1.  If  clearly  follows  fro-  the  construction  below  that 

* 

1  r  T  '  •  ,  :  '  =  0  ■For  u  \ 

X  t  f* -1  ^  1  «;n  ho 

* 

]  ir*  !j  t  ,  *  )  'r  '•)  % 

?rrr  c  .  '  ‘  ,:rcir';  r  't -■.*.•••  :  *.r>  ^'forcr  for  t  ho  caso  A  ~  1;  th^ 

~  1  /"> 

apr.f-ra  i  v '■  - •  t-r.  ♦-♦'.it  Sv  inq,  x  A  y  . 


-  >( 


For  some  w  £  f0,a/2],  lot. 


H(u'  =  .tf(ii),  u  .  r0,w) 

i 

\  y  f ( u ) ,  u  c  ( w, a-w> 
l 

f  {  u )  ,  u  c  (  a  -w ,  1  ) 

Pr,r  ...  ,  n  '1  H  ( r '  <*u  =  y  *a  fdu  -r  1  fdu  >  0,  since  y  <  v.  Therefore  there 

’  • 0  —  • 0  'a 

exists  a  value  c  (0,a/2)  for  which  /’  H(u)du  -  0  as  well.  For  w  =  w, 

let  “lx)  be  the  solution  of  (3.7)  satisfvino  U(0j  =  w  ^  .  I,  .it  (0,P  )  and 
(5*,  p*  +  T  )  he  the  intervals  on  which  F  t*  (w  ,  a  -  )  and  (a  -  w^,  a)  , 

respectively.  Then  ( 3.11)  for  P  =  r>* ,  T  —  T*  implies  (3.9a). 

:  s  f(u),  M  (9, a)  , 

Pef  i  r.  i  m  H  (  u)  = 

;  s^f(u>,  u  i  u,n  , 

we  g.-.e  that  (3.12)  implies  (3.6a);  and  if  "  is  the  solution  of  (3.P) 
satisfy.;'.';  !'(-X)  -  1,  (3.°h)  will  he  satisfied. 

r.*--r  p  --  n*  f  "  -  7 *  f  the  conclusion  follows  t yon  Theorem  3,  Now  let 

mi  "  to  anv  numbers  satisfvino  (3.13)  fX  still  eoual  to  1).  Then 
; „  .  y  •  It.  *  71  -  it*  «.  7* )  ,  v?e  see  that.  '3.11)  still  hold,  with  the 
.  ,  me  i  "r.e  4  t  Kore  shifted  to  the  riaht  by  the  amount  .  This 

.  ..  ...  :  .  1  .  \- 1  1  . *  d.-r  ;  -  ■  i  as  i  «-  tv  >  proof  of  oorol  1  ary  6.1. 

...  ;  1  '  !u)  1  i  ~f  f  I  .  Th.-»-  *  I'sl  inwirw  expressions  yi  eld 

’  ■  i 

r:  7  •  .  |  <  -  |  >  •  -  •  tv.,  .„a,,i  rr.~,.p»  r>*  enroll  ary  1: 


;•  7  - 


1  /2  “■  I 

B  =  (a  -  2w2  )  (F(  0  ,w^  )o  )  +  «,(F(a,D 

-1 

T  =  w  (F(a,1 ) ) 

•  ere  F(a,e)  2  (2  /^| f(u) !du) 1/2. 

Proof:  Let  c  <  0  be  such  that  U(c)  =  0.  Multiply  (3.7)  by  U'  an 

irtearate  from  c  to  0 ,  to  obtain 

.  1/2 

a’ (0)  >  p(0,w  )o 


e  ( 0  ,  B*  ) ,  0"  >  0,  so  U’(x)  >  D'lO),  and 


U  ( x)  >  w^  +  x  F(  0  ^  )0 


1/2 


Betti  no  x  =  B*,  U(B*)  =  a  -  ,  we  obtain 


1/2, 


<  (a  -  2w  KFfO.w^fl  ) 


(3.14) 


Working  from  the  other  end,  we 

*  *  * 

x  r  r  p  ,  B  +  T  1  , 


find  U'(B*  +  T  )  =  F(a,1),  and  for 


*  *  *  * 

U'(B  +  T  )  -  M(B  +T 


*  * 

-  x)  <  U' (x)  <  U' (P  +  T  1 


-vhere  -M  =  inf  f.  Integratinq  from  B  to 

*  *  * 

ry,ar.  ,,  (P  )  =  a  -  w  ,  n(B  +  T  )  =  a,  we  get 


p*  +  t*  and  using  the  facts 


-DP- 


*  1  *  * 

F( a, 1 )T  -  -  MT  <  w1  <  F(  a ,  1 )T 


hence 


-1  *  -1  2KW1  1/2 
w  F(a,1)  <  T  <  M  F  (  a  ,  1  )  [1  -  (1 - -)  ] 

(F(a,1))2 


<  2w  ( F( a , 1  )  ) 


(3.15) 


It  is  seen  from  (3.14)  and  (3.15)  that  the  values  of  B  and  T  given  in  the 


proposition  satisfy  (3.13), 


4 .  Properties  oe  solutions 


We  begin  with  a  monotonicity  theorem. 

Theorem  6.  Let  v‘  be  a  solution  of  Problem  I  in  which  s  ami  f 
satisfy,  respectively,  an!  Then  if  s  is  nondecreasing,  'f  is 

strictly  increasing. 

Proof .  Without  loss  of  aenerality  we  may  set  x  =  1.  Let 


x  =  infix  s  R  :  >  0  on  (x,“>)} 


and  let  us  assume  that  x_  >  Set 


v  2 


=  y{ x? ) .  Since  v'(x7)  =  0, 


y"(xo)  i>  0,  and  u  ^  a  is  a  solution  of  equation  (2.1),  it  follows  that 
v“2  £  (  0  ,  a)  and  (  x2  )  >  o.  Let 


x^  =  infix  <  x^  :  y'  <  0  on  (x,x2)} 


In  view  of  the  boundary  condition  at  x  =  x1  >  -  °°.  Set  ^  Xl ) 

Then  ~  ^  e  (a,1).  Denote  bv  x-|(y)  the  inverse  of  y(x)  on  (x1(X2>  and  by 
x2(y‘)  the  inverse  of  c(x)  on  (x2,<c). 

If  we  multiply  (2.1)  by  and  integrate  over  (x^,x2>  we  obtain 


,  s(x1  (  c)  )f  (y  )d  y  =  0  , 


(4.1  ) 


and  if  we  integrate  over 


I  x 


2 


,=°) : 


•1 


2 


s  I  x 


2 


n 


(4.2) 
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Equation  (4.1)  yields 


/*  s(Xl(v))  |f(v)  |dv* 


s(Xl  ( >^)  )  |  f  (  v)  I  d' 


and  hence,  by  the  monotonicity  of  s 


s(C ) 


f*  |f(v)|dv'  <  s(£)  /  1 1  f  ( |  dy? 

'  ^2  3 


whence  5  is  the  (unique)  zero  of  the  equation 


v  ( x )  -  a  =  0 


on  (  x-j ,  *2  )  •  Thus,  since  s  >  0 


J®  |f(  v')|d-'  <  J  1  t  f  (  V )  |  ds5  . 
"2  a 


Similarly,  (4.2)  yields 


,  "!■  s(x  (  v) )  I  f  (  v)  |  =  j1  s(  x  ( i^) )  |  f  (  <)  !  d^ 

v  2  ^  9  t 


an'i  ■  ipr.ctt 


s(r)  'a  |  f  (v  )  |  d  v  i  s(t)  /1|f(v')|dy- 
y'2  * 


i  r;  tho  (union**)  zoro  of  (4.3)  on 


'V1* 


Thus 


(4.5) 


fa„  |f<*)  Id*  >  J1  I  f  (sS>  |  dss  . 
J  *  '  a 

2 


Together  (4.4)  and  (4.5)  imply 


* 

J  1  |  f(*)  Id*  >  J1 1  f(*>  |d  *  , 

3  3 


which  would  mean  that  >  1.  This  is  impossible. 

Next,  we  turn  to  the  question  of  uniqueness,  we  shall  show  by  means  o c 
the  example  discussed  in  Problem  1,  that  we  generally  cannot  expect  a  unique 
solution  (see  also  [17]). 

Proposition  2.  Let  X  *  1  and  s  :  R  R+  be  defined  by 


r  S  x  <  0 

s(x)  a  J 

l  s2  X  >  0 

whence  0  <  s^  <  S2’  Then 

(i)  if  s^/52  =  Y*  Problem  I  has  precisely  one  solution 
(ii)  if  si/s2  <  Y»  Problem  I  has  precisely  two  solutions. 

Proof .  Let  v(x,a)  be  the  solution  of  equation  (2.1)  on  (-“,0)  which 
satisfies  the  boundary  conditions  v(-°=,a)  =  0,  v(0,a)  =  a.  Then  for  each 
»  c  [0,1],  v(*,a)  is  well  defined  and 


{v’(0,a)}2  =  -2Sl  J“  f  ( r)dr ,  n  4  a  4  1  . 

.similarly,  let  w(x,a)  be  the  solution  of  equation  (2.1)  which  satisfies  the 

boundary  conditions  v(0,ci)  =  a,  v(°°,a)  =  1.  This  solution  is  well  defined 

for  a  <  a  4  1 ,  whence  a  is  determined  by  f(r)dr  =  0,  and 

a 
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5 .  Bounded  donains. 

In  this  section  we  oonsider  the  problem 

r  u"  +  Xs(x)f(u)  =0  -1  <  x  <  1 

(III) 

V  u* (-1  )  *  u'  (  +  1  )  =  0  . 

Clearly  this  problem  has  three  trivial  solutions  u.|  =  0,  U2  =  a  and  =  1. 
u1  and  u^  are  stable  and  U2  is  unstable.  We  shall  inquire  into  the 
existence  of  stable  nontrivial  solutions. 

To  begin  with  we  shall  aeneralize  Theorem  4.  However,  before  we  can 
state  it  we  need  to  introduce  some  notation.  Consider  the  problem 

r  u"  +  x  f(u)  =  0 
u(0)  =  a,  u’ (1 )  =  0  . 

This  problem  has  a  strictly  increasing  solution  u(x,X)  for  X  >  X^  >  0.  Set 

...  ru  ( 1  , X )  ,  .  . 

w  (  X  )  =  f  (  r )  dr  X  >  X 

J  a  0 

Clearly  .j  (  X  1 '■-j  0  as  X  ^  X  and  ^(“l  =  {  f(r)dr 

U  3 

Next,  consider  the  problem 


u"  +  Xf(u)  =  0  (5.1) 

u ( 0 )  =  1 ,  u' (1)  =  0  .  (5.2) 


In  view  of  the  integral  condition  on  f  there  exists  a  X  >  0  such  that  if 
X  *  *  ,  this  prob1'>m  has  a  strictly  docreasinn  solution  (which  is  not 
unique ) . 
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T! 


Theorem  7.  Let  f  satisfy  and  let  s  be  of  the  form  (3.1)  in 

which  c  (0,1)  and  6  satisfies  assumptions  H  . 

Suppose  \  and  5  satisfy  the  following  inequalities 


XC2.o(X(1-rj2)  >  1  a2 

H 

X(1-£)2  >  max{X  ,X  } 


(5.3) 

(5.4) 


Then  there  exists  a  constant  y  such  that  if  y  >  y  ,  Problem  III  has  a 

* 

solution.  If  X  *  »,  then  u  •*  00 . 

Proof .  We  proceed  as  in  the  proof  of  Theorem  4  constructing  a  sub- 
solution  v  and  a  suoersnlution  w  such  that  v  <  w. 

ror  Vj  we  now  '■aXe  the  solutin  of  (2.1)  on  (^,1)  which  satisfies  the 
boundary  conditions  v  ^  ( \ )  =  a  and  ’ ( 1 )  =  0.  This  solution  exists  in  view 
of  (5.4)  and  it  can  be  seen  bv  a  scalina  araunent  that 


{ v  ’  f1:)}2  =  2  X  w  (  X  ( 1  ~C  )2 ) 


(5.5) 


be  nroof  of  Tbeoren  4  we  reauire  that 


v3  1  (\)  >  n/2K  f 


'r  wiM  (  )  yield*  (c.3  )  . 

;r  ;r*  i  x  r  r  c  + 


v  o  v ?  i  s  the  same  as  in  the  proof 


_  -)  r,  _ 


1 

i 

4 


The  supersolution  w  consists  again  of  two  nieces  and  w£. 

Again,  w2  “  ^  on  [“€#^1  but  is  now  the  solution  of  equation  (2.1)  on 

which  satisfies  the  boundary  conditions  (-1 )  =  0  and 

55  1  •  It  is  a  simple  matter  to  transform  this  problem  to  (5.1), 


6 .  Variable  migration  and/or  carrying  capacity 
We  now  consider  the  problem 


<k2 (x)V2 (x)u' ) 1  +  2 Xk2 ( x)V ( x) f ( u)  =  0,  x  e  R,  X  >  0 


u  (  _oo  )  =  0 ,  u  ( “ )  =  1  , 


where  f  satisfies  Hf,  k  and  V  are  in  C  and  k(x)  and  V(x) 

have  positive  lower  bounds.  As  we  shall  see,  many  of  the  results  about 
Problem  IV  can  be  deduced  from  the  results  obtained  about  Problem  I. 


,  .  def  rx  dr 

y  =  t<x)  =  /  — - j — 

k  (r)V  (r) 


Clearly  x  maps  R  onto  R  and  its  inverse  x  is  well  defined. 
Introducting  y  as  the  independent  variable  into  (6.1),  we  obtain 


u"  +  Xs(y)f(u)  =  0,  y  £  R 


where  the  primes  now  denote  dif ferentiation  with  respect  to  y,  and 


s(y)  =  2k4 ( x-1 (y) )V3(x-1 (y) ) 


*  4  3  *  4  3 

s1  =  inf  k  (x)V  (x)  and  s?  =  sup  k  (x)Y  (x) 

R  R 


Then  the  followinq  result  is  immediate  from  Theorem  2: 
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Theorem  8.  Let 


*  r1 

s  ;  f(u)du 
—  >  — 2 - 

s2  f(u)du 

Then  Problem  TV  has  no  solution. 

It  is  also  obvious  from  Proposition  1  that  Theorem  8  is  optimal.  Next, 
we  observe  that  for  any  x^,  x,  e  R. 

k4(x1)V3(x1)  -  k4(x2)V3(x2)  =  s(y1 )  -  s(y2)  , 

where  y^  =  x(x  J  (i  =  1,2).  Thus  if  k4(x)V3(x)  is  nonincreasing,  so  is 
s.  Whence  by  Theorem  2,  we  have  proved  the  following. 

Theorem  9.  Let  k4(x)V3(x)  be  a  nonincreasing  function  of  x.  Then 
Problem  TV  has  no  solution. 

Finally,  we  shall  show  that  like  s,  k^3  has  to  increase  somewhere 
sufficiently  fast  for  there  to  exist  a  solution.  Consider  the  problem 

r '22  2 

[k  (£x)V  (cx)u' ] '  +  2k  ( ex)V( ex) f (u)  =  0  , 

(IV  )  ' 

£ 

u  (-“>)=  0  ,  u  ( <*> )  =  1 

Then  writing  x'  =  ex,  we  can  reduce  this  problem  to  Problem  IV  with 
2 

\  =  1/e  ,  and  hence  to  Problem  T  with  y  =  t(x’),  and  s  defined  by 

(6.2)  •  Theorem  3  now  yields  the  desired  result. 

*  * 

Theorem  10.  There  exists  an  t  >  o  such  that  if  c  <  £  ,  Problem 
TV  has  no  solution. 

Next,  we  turn  to  t he  ones* ion  of  existence.  It  is  here  that  the  results 
for  Problem  TV  heenme  somewhat  different,  and  in  fact  stronger. 
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To  hecjin  with,  we  introduce  the  two-parameter  family  of  functions 


V(x,u,£ )  =  {1  +  u6(x,\)}  ,  u  >  0,  £  >  0  , 


where  k  £  C  (R)  is  as  in  section  3. 

Theorem  11^,.  Let  f  satisfy  Hf,  and  let  V  be  of  the  form  (6.3)  in 
which  satisfies  assumptions  H ^ .  Let  k  e  C^(R)  satisfy  0  <  k(x)  <  1, 

and  k(x)  3  1  for  |x|  >  .  Then  for  each  >.  >  0 ,  £  >  0 ,  there  exists  a 

★  * 

li  >0  such  that  if  u  >  u  ,  Problem  IV  possesses  a  cline. 

Theorem  11,..  In  the  hyootheses  of  Theorem  1 1y  and  in  (6.3),  replace  the 
symbol  v  by  k  and  k  by  V.  Then  under  the  new  hypotheses.  Problem  IV 
possesses  a  cline. 

Proof .  We  aive  the  proof  for  Theorem  1 1V  only,  as  that  for  Theorem  11^ 


is  the  same.  Throuqhout,  we  shall  keep  £  fixed,  and  we  shall  write 
V  =  V(x,;j)  and  0  =  (x).  Define 


r'+(u)  =  -0  - -  *  4  vTrTuT  =  5  +  U  /5  6(r,dr  ; 

k  (  r )V  ( r  ,  u) 


(U)  = 


dr  _  ,0  dr  .  ,0 

— - >  ,-r  7(777T  =  *  +  u  /_e  e(r)dr 

k  (r)V (r,y)  ’  ,W  ^ 


f-)  ♦  ™  as 


•  rv-p  proof  of  Theorem  4,  we  construct  a  subsolution  v  which 
■f  three  Pie,-,.-,:  ,  vy,  and  v^.  The  last  piece,  V3 ,  is  the 


■'re,  wit' 


replaced  bv 


Y  r.  <t;r)  t". 


-  p  1  2  fl  -1 

1  )  >  —  a  (  ^  f ( r ) dr ) 


I 


where  n*  =  ni(u.]>.  Then  for  u  >  y  , 

1  +  +  — 

v3<n  )  >  a/(n  +  n  ) 


We  now  define  vi<y)  =  0  for  -“  <  y  <  p.  The  composite  function  v  is  the 
desired  subsolution. 

For  the  supersolution  w,  we  can  take  the  same  function  as  in  the  proof 
of  Theorem  4,  with  -£  replaced  by  -ri~(y). 

Thus  if  u  >  u  =  maxiul»U2//  there  exists  a  subsolution  v  and  a 
supersolution  w  such  that  v  <  w  on  R.  This  implies  the  existence  of  a 
stable  cline. 

Remarks .  1.  In  the  case  k  =  1,  it  follows  from  Theorem  P  that  the 

parameter  u  in  Theorem  1 1V  must  satisfy 


♦ 
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7 .  It  is  clear  from  (6.4)  that  if  A  +■  0,  then  y  03  and  hence 

u  ♦ 

Next,  we  <:  sider  (6.1)  on  a  hounded  domain: 

f (k2 (x)V2 ( x)u' ) '  +  Ak2(x)V(x)f (u)  “0,  -1  <  x  <  1  , 

(V)  < 

{ u 1 (-1 )  =  0,  u' (  +  1 )  »  0 

Theorem  12.  Let  f  satisfy  Hf  and  let  V  he  of  the  form  (6.3)  with 
0  <  5  <  1  in  which  6  satisfies  assumptions  Ha.  Let  k  e  C1)*)  satisfv 

D 

0  <  k(x)  <  1,  and  k(x)  s  1  for  |x)  >  £.  Then  provided 

A(1-S)2  >  max{AQ,A  } 

*  * 

there  exists  a  u  >  0  such  that  if  y  >  u  .  Problem  V  has  a  cline. 

* 

Recall  that  and  A  were  defined  in  section  5. 

Proof .  Following  the  proof  of  Theorem  11,  we  transform  equation  (6.1) 
on  (-1.1)  to  equation  (2.1)  on  the  interval  I  =  I1  u  I2  u  Ij,  where 

i1  =  (-n  -  (1-0,  -  n  ),  i2  *  C-n  ,n+l,  i3  =  (n  ,n+  +  d-O) 

On  Ii  and  I3,  we  choose  the  same  super,  resp  subsolution  as  in  the  proof 
of  Theorem  7,  and  on  Ij  we  choose  them  as  in  the  proof  of  Theorem  11. 
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7 .  Clines  on  a  bounded  habitat  II. 

In  this  section  we  prove  the  existence  of  dines  in  bounded  habitats  by  a 
method  which  is  different  from  the  one  used  in  the  previous  sections.  There 
our  main  tool  was  the  maximum  principle;  here  it  is  a  variational  argument  due 
to  Matano  [13].  He  used  it  to  prove  that  the  problem 

f  Au  +  f  (u)  =  0  in  si 
l  ^  =  °  on  3^  , 

where  fi  is  a  bounded  domain  in  rNn  >  1)  with  smooth  boundary  3fi,  and 
3/3n  is  the  outward  normal  derivation  on  SQ,  could  have  dines  if 
consisted  of  subdomains  £  which  were  connected  by  sufficiently  narrow 
passages.  The  idea  was  that  in  this  manner,  the  passages  offered  an  ob¬ 
struction  to  the  migration  of  the  individuals  between  the  subdomains,  thus 
allowing  them  to  maintain  different  states. 

The  problem  considered  in  section  6  is  analogous  to  the  one  considered  by 
Matano:  the  obstruction  now  being  caused  by  a  diminished  mobility  V  or 

carrying  capacity  k. 

V.'e  consider  the  aeneral  problem 

(Dfx)u')’  +  Xs(x)f(u)  =0  X  E  fl  =  (-1,1) 

(VI)  (7.1) 

u’ (-1  )  =  0,  u' (  +  1  )  =  0 

and  we  assume  that  the  function  f  satisfies,  in  addition  to  the  hypotheses 
!ir  the  normalisation: 

*  ? 

H ,  fr.-a)f(u)  (u-a)  for  0  <  u  <  1 
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About  D  and  s  we  make  the  assumptions: 


1+y  — 

.  P,s  s  C  (fi)  for  some  v  e  (0,1] 


A0.  0 

<  min{D(x)  : 

X  Z  il} 

and 

D(x) 

<  1 

on 

n 

0 

<  min  is ( x)  : 

x  e  P.) 

and 

s(  x) 

<  1 

on 

a  , 

A3.  1? 

(x)  =  1  and 

s(x)  =  1 

on 

°2 

t 

==  (-1,- 

- , )  and  1 

=  <S,1> 

for 

some 

5  £ 

(0,1) 

.  Still  following 

Matano,  we  define 

R  ( - ,  +  ]  =  {.;  £  C2(u)  C1  (I!)  :  /  (c-a)dx  <  0,  /  (;-a)dx  >  0} 

1  “2 

and  we  introduce  the  functional  J  :  H'(fi)  +  R: 

J(C)  =  fr[~  D<x){;' (x)}2  -  Xs(x)F(C))]dx  , 

where 


F(U  =  /C  f  ( r)dr  . 

3 


Recall  that 


/ 1 f ( r ) dr 
v  a _ 

-jrpf(r)dr 


F  (  1  ) 
F(0) 


Lomma .  Let  f  satisfy  the  hypotheses  Hf  ana  Hf*,  and  let  D  and 
s  satisfy  the  assumotins  A1  -  3.  Suppose  there  exists  a  function 
w  t  •!*(..)  such  that 
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2 

J(w)  <  XF(0 ) { ( 1 -£ )min( 1 ,  — — — )  -  2}  . 

(1-U  X 


Then  Problem  VI  has  a  stable  solution  in  the  set  R[-,+].  The  proof  of  this 
result  is  nearly  identical  to  that  of  Theorem  6.2  in  [13],  whence  we  omit  it. 

To  obtain  explicit  conditions  on  £ ,  X ,  D  and  s,  which  ensure  that 
(7.2)  is  satisfied,  we  consider  the  function 


We  find  that  J(w) 


w(  x)  = 


satisfies  (7.2)  if  £, 


-1  <  x  <  -5 

<  <  x  <  C 

5  <  x  <  1  . 

X  and  D  satisfy  the  inequality 


~  D ( x ) dx  <  41gF(0)[(1-^)min(-  1 )  -  1  +  y  -  d+y)?l  •  (7.3) 

S  (1-5)  X 

Thus  we  have  proved  the  following  existence  theorem. 

Theorem  13.  Let  f  satisfy  hypotheses  and  and  D  and  s 

assumptions  A1  -  3.  Then,  if  X,  £  and  D  satisfy  the  inequality  (7.3), 
Problem  VI  has  a  cline. 

Remarks .  1.  Note  that  beyond  the  assumptins  A1  -  3,  no  conditions  on 

s  are  required. 

2.  For  the  equation  discussed  in  section  6,  the  condition  for  the 
existence  of  a  cline  becomes  an  upper  bound  for  the  integral 


which  demonstrates  again  that  in  this  context,  the  roles  of  V  and  k  are 


interchangeable. 


s<  IjPt  7=1,  ana  C  -  1/4.  Then  (7.3)  becomes 


?  ty*  D(x)dx  <  J  F(O)RtA)  0  <  X  <  2  X0 


whence  Xn  =  <4ir/3)  and 


R(X)  = 


j'  X  for  0  <  X  <  XQ 

'i  3X  -  2X  for  x0  <  x  <  |  X0 
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